Draft version September 16, 2011 

Preprint typeset using I^'T^]X style eniulatcapj v. 11/10/09 



O 
(N 

CD 
00 



o 

I 

o 



(N 
> 

in 
in 
vn 

m 

i> 

o 



% 



VELOCITY EVOLUTION AND THE INTRINSIC COLOR OF TYPE lA SUPERNOVAE 
Ryan J. Foley^'^, Nathan E. Sanders\ and Robert P. Kirshner^ 

Draft version September 16, 2011 

ABSTRACT 

To understand how best to use observations of Type la supernovae (SNe la) to obtain precise 
and accurate distances, we investigate the relations between spectra of SNe la and their intrinsic 
colors. Using a sample of 1630 optical spectra of 255 SNe, based primarily on data from the CfA 
Supernova Program, we examine how the velocity evolution and line strengths of Si II A6355 and 
Ca II H&K are related to the B — V color at peak brightness. We find that the maximum-light 
velocity of Si II A6355 and Ca II H&K and the maximum-light pseudo-equivalent width of Si II A6355 
are correlated with intrinsic color, with intrinsic color having a linear relation with the Si II A6355 
measurements. Ca II H&K does not have a linear relation with intrinsic color, but lower-velocity SNe 
tend to be intrinsically bluer. Combining the spectroscopic measurements does not improve intrinsic 
color inference. The intrinsic color scatter is larger for higher- velocity SNe la — even after removing 
a linear trend with velocity — indicating that lower-velocity SNe la are more "standard crayons." 
Employing information derived from SN la spectra has the potential to improve the measurements of 
extragalactic distances and the cosmological properties inferred from them. 
Subject headings: supernovae; general — distance scale — dust, extinction 



1. INTRODUCTION 

Type la supernovae (SNe la) are very good distance 
indicators after making an empirical correction based 
on their light-cu rve shape and color (jPhillipsI 119931 : 
iRiess et al.l 119961 ) . Using these relations, large SN la 
samples h ave a precision o f ^8% for distance estimates 
(e.g., jHicken et al.l l2009aD . This level of precision is 
adequate to determi ne that the expa. i ision of the Uni- 
verse is accelerating (jRiess et al.lll998l: iPerlmutter et al.l 
Il999| ). and constrain t he equatio n-of-state parame ter of 
dark energy (IWood-Vas ev et al.l 2007; Riess et al.l 



Hicken et al. 2009b; Kessler et al. 2009; AmanuUa h et al.l 
2010: Conlev et al. 2011 ; Sullivan et al. 20 11). SNe la are 
even better dis tance indicators in the r est-frame near- 
infrared (NIR) (jMandel et al.|[200l[20lTI ): however, low- 
redshift NIR samples are small and high-redshift NIR 
samples do not yet exist^ 

Decades ago. Branch (Il987| ) showed that not all SNe la 
have the same ejecta velocity, as pr obed by the blueshift 
of spectral features. [Branch et al.l ([1988,) noted that a 
sample of SNe la had a very broad distribution of Si II 
A6355 velocity near maximum brightness and concluded; 
"SNe la are not observationally homogeneous. The only 
way to maintain that they are physically homogeneous 
would be to postulate that the y are identical but asym- 
metrical." iBenetti et al.l (|2005f) showed that the velocity 
gradient for the Si II A6355 feature did not correlate with 
light-curve shape for spectroscopically normal SNe la. 
This observation provided further proof that SNe la were 
heterogeneous, but also showed that some properties of 
a SN la did not depend on the width of its light curve. 

Several studies have shown that SN la spectral prop- 
erties are related t o their photometric properties (e.g., 
iNugent et al.l 119951 ). Some studies examined the corre- 
lation between Hubble residuals after light-curve shape 
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correction and spectros copic properties to produce a 
smaller Hubble scatter (Folcv et al. 2008b; 'Bailev et al.l 
i2009; Blondin et al. 2011b; ?; Nordin et al. 2011b). How- 
ever, using this simple approach with only optical spec- 
tra provides only mh ior gains over photometry alone 
(IBlondin et al.ll2011bl) . 

In a previous study, iFolev & KasenI (|201lL hereafter 
FKll) showed that the maximum-light intrinsic color of 
SNe la is highly correlated with their ejecta velocity as 
probed by the Si II A6355 feature, usi ii. Accounting for 
this correlation both substantially improves the preci- 
sion and reduces potential biases of SN la distance mea- 
surements. Theoretically, this is understood as increased 
line blanketing in the higher- velocity SNe, depressing the 
_B-band flux relative to the T^-band flux and causing a 
redder B — V color. Another theoretical expectation is 
that higher-velocity SNe la have more absorption from 
lower excitation (e.g., Fe II vs. Fe III) lines, which are 
distributed such that higher-velocity SNe la should have 
more absorption in the near-ultraviolet. Nonetheless, not 
all SN la model s show the faste r-redder relation found in 
the data ( Blondin et al.ll2011aD . 

Although not necessarily the cause of the relation be- 
tween intrinsic color and vgi ii, the observatio ns are nat- 
urally explained by as ymmetric exp l osions (jFKllt ). as 
originally suggested by iBranch et al.l (|1988f ) . An asym- 
metric explosion model was invoked to explain a strik- 
ing relation between the velocity gradient of the Si II 
A6355 feature, {jsi n , and the velocity of fset of nebular 
lines a t late times (jMaeda et al.ll2010l ). iLeonard et al.l 
(|2005f ) found that SNe la with high- velocity features had 
stronger Si II A63 55 line polarization than those with- 
out the features. iMaund et al.l ()2010( ) expanded upon 
this work and found a correlation between the amount 
of polarization in the Si II feature and vsi n, suggest- 
ing that wsi II is an excellent probe of asymmetry of the 
outer layers of t he SN ejecta. Since vsi n is an excellent 
proxy for wsi u (jWang et al.ll2009al . hereafter W09), one 
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expects fsi II, wsi ii, intrinsic color at maximum bright- 
ness, nebular line velocity offse ts, and Si II polariz ation 
should all be highly correlated. iMaeda et alj (|201lD con- 
firmed one of these assumptions, showing that nebular 
line velocity offsets are correlated with maximum-light 
color. Of all of these observables, wsi ii is the easiest to 
obtain at low redshift (and is always obtained if iisi ii 
or polarization measurements are made) and is the only 
observable known to correlate with intrinsic color that 
we can realistically obtain for high- 2; SNe. 

Besides line velocity, which is determined by the min- 
imum of an absorption feature (corresponding to its 
maximum absorption), one might expect that the total 
amount of absorption, as measured by either the full- 
width at half maximum (FWHM) or pseudo-equivalent 
width (pEW) of a feature, to strongly correlate with in- 
trinsic color. Both the FWHM and pEW are alternative 
kinematic probes, measuring the span of the absorbing 
material in velocity space and may have more physical 
motivation for a correlation with intrinsic color than ab- 
sorption velocity. 

Although Si II A6355 is a clean, isolated, strong feature 
in SN la spectra, there is no guarantee that its properties 
are the most informative for determining intrinsic color. 
Furthermore, Si II A6355 redshifts out of the optical win- 
dow at 2: « 0.4. Therefore, any relation determined with 
Si II A6355 cannot be applied to most high-z SNe la. 
Another strong feature in SN la spectra that has a large 
velocity distribution and st rong velocity evolution with 
time is Ca II H&K. iFKllI suggested that this feature, 
which is much bluer than Si II A6355 and can be ob- 
served in the optical window up to z « 1.2 , may be an 
alternative way to determine intrinsic color. iFKllI found 
that the model spectra of IKasen fc Plewal (|2007l ) show 
a strong relation between intrinsic color and fca h&k, 
similar to what was found with vsi u- 

Current and previous high-2 SN la samples typi- 
cally have a single near-maximum b rightness observer- 
frame opti cal spectrum per SN (e.g. , [Zheng et al]|2008l : 
iFolev et al . 2009; Walk er et al.ll20lil ). making measure- 
ments of V impossible. Many current and future surveys 
intend to use photometric-only samples of SNe for cos- 
mological studies where only a small subsample of SNe 
will have spectra. These samples will also not have the 
spectral series necessary to determine velocity gradients 
for high-2; SNe la. Moreover, if intrinsic color can be 
determined by a single spectrum, telescope time can be 
used to make this measurement for more SNe than one 
could if a time series is requi red. 

Since IFK 111 used the lW09l sample, which did not pro- 
vide spectra or velocity measurements, but only a des- 
ignation of "Normal" or "High- Velocity" for SNe, FKllI 
was unable to investigate exactly h ow in trinsic color de- 
pends on ejecta velocity. Similarly, IW09I did not provide 
any FWHM or pEW measurements or any measurements 
for Ca II H&K. In this paper, we use the large and homo- 
geneous CfA sample of SN la spectra (Blondin et al., in 
prep.), supplemented by literature data (Section [2|), to 
construct a homogeneous sample which reduces poten- 
tial systematic differences, determine a relation between 
V and v{t) for the two strongest features in a maximum- 
hght SN la spectrum, Si II A6355 and Ca II H&K (Sec- 
tion [3]) , use the same data to determine a relation be- 
tween pEW for these features and their gradients (Sec- 



tion 2]), use those relations to measure maximum- light 
values for a large sample of SNe, allowing a direct com- 
parison of these SNe (Sections |3] and S]) , and reexamine 
the relation between ejecta kinematics — as determined 
by both the width and velocity of the minimum of the ab- 
sorption for Si II A6355 and Ca II H&K — and intrinsic 
color (Section[5]). We find that the maximum-light veloc- 
ity of Si II A6355 and Ca II H&K and the maximum-light 
pEW of Si II A6355 are correlated with intrinsic color. 
The intrinsic color scatter is larger for higher-velocity 
SNe la — even after removing a linear trend with ve- 
locity — indicating that lower- velocity SNe la are more 
"standard crayons. "0 We discuss these results and con- 
clude in Section [6l 

2. DATA 

For our study of how SN la intrinsic color is related to 
kinematics, we require a large sample of SNe la with both 
photometry (to measure colors) and spectra (to measure 
kinematic probes). Below we describe the data used for 
our sample. 

Ad ditionally, we use values reported by IW091 Since 
IW09I does not report velocity measurements (only a 
designation of "Normal" or "High- Velocity" ) , it is not 
sufficient for use as a spectroscopic sample. They 
do report some relevant photometric measurements 
(specifically My at peak, Ami5{B), and i?max — 
t4iax)- How ever, the photometric data come from var- 
ied sources ("Hamuv ct al.l Il996al iHamuv et a l.' '1996bt 
CfAl; Ricss ct al. 1999 ; al l compi led by Rein dl et aT] 
l2005t CfA2: lJhaet all 12006: CfA3; Hick en et al.ll2009ar 
and preliminary LOSS ligh t curves later published by 
iGaneshalingam et al.l I2010D . and there are potentially 
systematic differences between each data set. A com- 
pilation of several photometric sources can have major 
differences, even after placing all photometry on a stan- 
dard system. At the very least, S corrections are typ- 
ically nficessar5;_to_a£hieve agreement at the 0.01 mag 
level (jStritzinger et al.ll2002D . 

Although we use the data presented by IW09I for an 
initial analysis (see Section [5]) , we do not combine data 
from , W09i with other data for our final analysis. 

2.1. Photometry 

Our photometric d ata comes fro m two sources: CfA3 
(iHicken et al.l l2009aD and LOSS (jGaneshalingam et al.l 
I2OIOD . Each sample is large, and there is substantial 
overlap in SNe between the two groups. Each group 
derived parameters from their light curves. Specifi- 
cally, they each measure i 3max, Vmax, and Ami5(_B). 
IGaneshalingam et al.l ()2010D used a combination of tem- 
plate fitting (when the data are good and similar to nor- 
mal SN la light curves) and polynomial fitting (when 
template fitting fails). Hicken ct al. (2009a) only used 
polynomial fitting. The template fitting is more sophis- 
ticated and should produce better measurements. We 
therefore use LOSS-derived parameters for SNe la that 
have both LOSS and CfA3 data. 

For 40, 42, and 29 SNe la in common for CfA^ and 



LOSS, we find average offsets of {B] 



CfA 
max 



B, 



LOSS 



^ See IPeek &: Graves! 1)20101 ) for a description of this nomencla- 
ture. 

■* We correct the CfA3 SN 2008bf photometry for an incorrect 
calibration which affected the light curves by 0.07 mag in B and 
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0.020±0.010mag, {V, 



CfA 



K 



LOSS 



) = -0.006±0.008mag, 



and (Amfl^ - l\m\^^'^) = 0.008 ± 0.017 mag, respec- 
tively. Furthermore, for 40 SNe in common, we find 

((Bmax - K.ax)^f^ - (Sn,ax - K.ax )^°^^) = 0.025 ± 

0.011 mag. Considering that iFKllI found an offset in 
Smax— Vmax of ~0.1 mag for low and high- velocity SNe la, 
this offset could cause enough scatter to reduce the sig- 
nificance of potential correlati ons. 

Gancsh alingam et al.l (j2010( ) examined the light curves 
in common between the LOSS and CfA3 sample. In- 
terpolating LOSS light curves, they examined differ- 
ences between individual light-curve points. They found 
average offsets of 0.011 ± 0.006 (0.016 ± 0.005) and 
-0.006 ± 0.004 (-0.010 ± 0.004) mag for B and V us- 
ing only data brighter than magnitude 18 (all data), re- 
spectively. These measurements are consistent with our 
finding for the derived parameters. Confirmation with 
the authors indicates that extinction corrections and K 
corrections were performed in the same way (M. Gane- 
shalingam & M. Hicken, private comm.), leaving only S 
corrections as a possible difference. Since a SN spectrum 
changes with time, the S corrections should also change 
with time. Therefore, making corrections to derived 
peak magnitudes using the offset w e found above is likely 
more accurate than that found by iGaneshalingam et al.l 
()2010D . which used the full light curve. 

Considering the relatively large offsets, we correct the 
CfA3 photometry to match the LOSS photometry, which 
is already the default photometry when there are data 
from both samples. 

A ll derived values presented b y iHicken et all ()2009aD 
and IGaneshalingam et al.l ()2010[ ) are both K-correct and 
deredden for Milky Way reddening. We present the com- 
bined derived light-curve parameters in Table [T] 



2.2. Spectroscopy 
2.2.1. CfA Spectral Sample 

Over the last two decades, the CfA Supernova Pro- 
gram has observed him dreds of SNe la, mo stly with 
the FAST spectrograph (jFabricant et al.l ll998') mounted 
on the 1.5 m telescope at the F. L. Whipple Observa- 
tory . The data have been reduced in a consis tent man- 
ner (jMatheson et al.ll2008l: [Blondin et al.ll2011b: Blondin 
et al., in prep.), producing well-calibrated spectra. 

For SNe la in the sample with a measured time of max- 
imum brightness from light curves, vsi ii and wca h&k 
have been measured (Blondin et al., in prep.). Briefiy, 
this is achieved by first generating a smoo thed spectrum 
using an inverse-variance Gaussian filter (jBlondin et al.l 
l2006f) . and the wavelength of maximum absorption in 
the sm oothed spectrum is u sed to determine the veloc- 
ity (see lBlondin et aDl2011bl for details). When possible 
the spectra are used to me asure the FWHM and pEW 
for the features as well f see lBlondin et aLll2011bl for de- 
tails). The measurements for each spectrum are listed 
in Table [21 The reference for each spectrum is listed in 
Table [21 but measurements in all cases were obtained by 
Blondin et al., in prep. 

0.13 mag in V . After correction, the peak magnitudes are Bmax = 
15.65 ±0.04 mag and Vmax = 15.81 ±0.04 mag (M. Hicken, private 
comm.). 



2.2.2. Literature Spectral Sample 

Many SN la spectra obtained by various groups us- 
ing various telescopes, instruments, and reduction meth- 
ods have been published over time. Although the fi- 
delity of the spectra may not be high or uniform as 
for spectra from the CfA sample, these data are use- 
ful for expanding our sample. When publicly available, 
usually through the SUSPECT database^, and if the 
spectra cover the appropriate wavelength ranges, wsi ii, 
pEW(Si II), VQa. H&K, and pEW(Ca II) have been mea- 
sured for these SNe. Since the flux calibration is not crit- 
ical for the velocity measurements, we expect the largest 
errors to be from poor wavelength calibration; however, 
the accuracy necessary for our measurements is easily 
obtained even if mistakes are made while observing or 
in the reduction process. Comparison of multiple spec- 
tra from the same SN both within a dataset and with 
the CfA sample show that the measurements are con- 
sistent within the measurement uncertainty (typically 
100 kms^^). The pEW measurements are more affected 
by poor flux calibration. We have examined the individ- 
ual spectra and have rejected spectra where we believe 
that the flux calibration is poor. The measurements for 
each spectrum are listed in Table [21 

Between the CfA sample and Literature sample, there 
are 1630 usi ii, 1630 pEW(Si II), 1192 wca h&k, and 1234 
pEW(Ca II) measurements for 255, 255, 192, and 211 
SNe la, respectively. For the majority of this study, we 
will restrict our sample to SNe with 1 < Amj5(i?) < 
1.5 mag to match the final sample criteria of lFKllL With 
this restriction, there are 939, 939, 685, and 708 measure- 
ments for 141, 141, 109, and 119 SNe, respectively. 

2.2.3. Telegraphic Spectral Sample 

Despite the CfA sample being quit e larg e, there are 
several SNe in the CfA3, LOSS, and [W09| photometry 
samples that do not have a CfA or publicly available 
spectrum (or one with —6 < i < 10 d). However, every 
SN in these samples was spectroscopically classified in 
the lAU Circulars or CBETs. In these reports, the date 
of observation is always reported, and for most SNe la 
classified near maximum brightness, fsi ii is reported. 
Although these measurements may not be of the same 
quality — and are certainly not systematically measured 
in the same way, as we have done with the CfA and 
public spectra — the uncertainty in the velocity mea- 
surement should be relatively small. These additional 
measurements from the Circulars and Telegrams can sig- 
nificantly expand our final sample. 

For most SNe, we are able to determine the phase of 
each spectrum using the time of maximum reported by 
IHicken et all (|2009b)Q and Gancshalingam et al. (2010). 
We present the data for these SNe/spectra in Table [3l 

2.2.4. Fll Sample 

The intersection of the CfA, Literature, and Tele- 
graphic spectral samples with the LOSS and CfA3 pho- 
tometric samples is the sample used for the analysis here. 



^ http://suspect.nhn. ou.edu/'suspect/ 
^ These values are 



reported 



at 



|http : //mm . cf a . harvard . edu/supernova/Cf A3/sn . tBmax . mlcsl7 ■ txt] 
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TABLE 1 
Object Information 



SN 



2CMB 



(mag) 



(mag) 



Bn 



lax V'rn 

(mag) 



"Si II 
(10=* kms-i) 



pEWo(Sill) «g,H&K' pEWo(Call) 

(A) (103 kms-i) (A) 



IQSIB'^ 


0.0036 


11.90 


(0.05) 


1.08 


0.03) 


0.01 (0.06) 


1986G 


0.0027 


11.12 


(0.06) 


1.75 


0.00) 


0.88 (0.09) 


1989B 


0.0035 


11.89 


(0.05) 


1.16 


0.02) 


0.32 (0.10) 


1990N<i 


0.0051 


12.66 


(0.04) 


1.11 


0.01) 


0.01 (0.06) 


1990O 


0.0306 






0.95 


0.02) 




1991M 


0.0076 












1991T 


0.0069 


11.46 


(0.02) 


0.93 


0.01) 


0.12 (0.03) 


1991bg 


0.0045 


13.86 


(0.04) 


1.94 


0.00) 


0.71 (0.06) 


1992A 


0.0058 


12.49 


(0.01) 


1.48 


0.00) 


0.02 (0.02) 


1992ag 


0.0259 


16.17 


(0.06) 


1.11 


0.09) 


0.01 (0.07) 


1992al 


0.0141 


14.56 


(0.04) 


1.10 


0.08) 


-0.13 (0.04) 


1993ac 


0.0493 






1.19 


0.10) 




1994D 


0.0041 


11.83 


(0.02) 


1.42 


0.00) 


-0.04 (0.04) 


1994M 


0.0243 






1.35 


0.03) 




1994S 


0.0160 


14.79 


(0.06) 


1.02 


0.00) 


-0.02 (0.08) 


1994T 


0.0356 


17.15 


(0.04) 


1.51 


0.08) 


0.18 (0.05) 


1994ae'^ 


0.0062 


13.00 


(0.03) 


1.09 


0.01) 


-0.04 (0.07) 


1995D 


0.0078 


13.26 


(0.05) 


0.98 


0.01) 


-0.10 (0.07) 


1995E 


0.0121 


15.98 


(0.05) 


1.18 


0.01) 


0.70 (0.07) 


1995ald 


0.0067 


13.20 


(0.05) 


0.93 


0.03) 


0.13 (0.07) 


1996C 


0.0275 






0.97 


0.01) 




1996X 


0.0078 


13.03 


(0.04) 


1.24 


0.00) 


-0.05 (0.05) 


1996Z 


0.0085 






1.22 


0.10) 




1996ai 


0.0037 






0.96 


0.02) 




1996bk 


0.0070 






1.78 


0.01) 




1996bl 


0.0348 






1.10 


0.08) 




1996bo 


0.0163 


15.51 


(0.04) 


1.20 


0.01) 


0.31 (0.06) 


1997E 


0.0133 


15.08 


(0.07) 


1.46 


0.02) 


0.03 (0.09) 


1997Y 


0.0165 






1.17 


0.02) 




1997bp 


0.0094 


13.74 


(0.03) 


0.96 


0.03) 


0.10 (0.04) 


1997bq 


0.0095 






1.01 


0.05) 




1997br 


0.0081 


13.43 


(0.08) 


1.10 


0.03) 


0.16 (0.13) 


1997cn 


0.0176 






1.90 


0.05) 




1997do 


0.0105 






0.94 


0.04) 




1997dt 


0.0061 






1.04 


0.15) 





-11.87 
-10.48 
-10.56 
-9.38 
-11.99 
-12.81 
-9.62 
-10.04 
-14.04 
-11.74 
-11.26 
-13.21 
-11.24 
-12.36 
-10.72 
-13.43 
-10.98 
-10.33 
-11.16 
-13.01 
-10.69 
-11.17 
-12.14 
-10.82 
-12.16 
-12.34 
-12.25 
-12.01 
-11.03 
-15.62 
-14.24 
-11.60 
-9.61 
-13.38 
-11.33 



(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.22) 
(0.30) 
(0.30) 
(0.22) 
(0.24) 
(0.27) 
(0.24) 
(0.30) 
(0.23) 
(0.23) 
(0.24) 
(0.23) 
(0.22) 
(0.22) 
(0.23) 
(0.22) 
(0.22) 
(0.24) 
(0.24) 
(0.23) 
(0.22) 
(0.22) 
(0.23) 
(0.22) 
(0.22) 
(0.25) 
(0.22) 



128.8 (6.0) 

121.1 (6.0) 
84.7 (6.0) 



-15.31 (0.57) 112.7 (6.9) 
-11.65 (0.57) 106.4 (6.9) 



109.6 (6.0) -17.22 (0.57) 98.9 (6.9) 



144.6 (6.0) 
100.6 (6.0) 
118.6 (6.1) 
85.0 (6.0) 

81.6 (6.0) 

98.9 (6.0) 



90.1 (6.0) 
110.6 (6.0) 



89.0 (6.1) 

133.5 (6.1) 

120.6 (6.0) 
101.1 (6.0) 

156.0 (6.0) 



82.6 (6.0) 



-17.93 (0.57) 
-11.53 (0.58) 



-12.67 (0.59) 
-15.79 (1.08) 



-11.70 (0.57) 
-16.53 (0.57) 



103.0 (6.9) 
92.4 (6.9) 

117.8 (6.9) 

110.0 (6.8) 

103.0 (7.8) 



80.7 (6.9) 
61.6 (6.9) 



-15.53 



(0.57) 51.8 



(6.9) 



Note. — As described in the text, all photometry has been places on the LOSS system. Uncertainties are in parantheses. 
^ C orrected for Milky Way extin ction using the reddening maps of lSchlegel et all II1998I ). as performed bv lHicken et al.l Il2009al 'l 
and IGaneshalingam et al.l | |201(]| ) . 

When color inforamation is available, Amis (B) has been corrected for host-galaxy extinction usin g the method similar to that 
of Phillips ct al. (1999); Ami5(B) = Ami5_ obs(-B) + 0.1 x ((Bmax -Vmax) -1-0.081) (see Section [52] The values for Bmax - Knax 
include a correction for Milky Way extinction, so that correction is inherent in this process. 
'^ Values given for the "red" component. 

"* Cepheid object with distance modulus from lRiess et afl 1120111 ). The redshifts listed for these SNe yield appropriate distance 
moduli using Hq = 70.5 km s~^ Mpc"'^. 



Similar to lFKll | . we exclude SN 2006bt because of its pe- 
culiar nature (JFolev et "aL ll2010[). We also exclude from 
our final analysis SNe with E{B — V)mw > 0.5 mag, cor- 
responding to SNe 1999ek {E{B - y)Mw = 0. 561 m ag) 
and 20 061f {E{B - 1/)mw = 0.954 mag). Both |WOi and 
IFKIII excluded SN 20061f from their analyses. We call 
the final sample of SNe with photometric and spectro- 
scopic measurements the "Fll" sample. 

3. VELOCITY GRADIENTS & MAXIMUM-LIGHT EJECTA 
VELOCITY 

The temporal coverage of the Fll spectra (in partic- 
ular, the CfA data) provide an ideal sample for mea- 
suring velocity evolution with time. Previous stud- 
ies have examined relatively small samples of nearby 
SNe la and their sp ectra to determine velocity gradients 
for th ose SNe (e.g. JBenetti et all 120051 : iHachinger et alJ 
I2006D . Additionally, ejecta velociti es have been mea- 
sured for some h i gh-z SNe la (iHook et al.l 120051: 
Blondin et al.' '2006': 'Brondcr et al.' '20081: iNordin et al.l 
2011a: iKonishi ct al.. i2011.) . WO 9. separated their sam- 



ple into two subsamples based on ejecta velocity near 
maximum brightness. Since SNe la with high velocities 
near maximum brightness tend to have high velocity gra- 
dients, this separation is similar to separating based on 
velocity gradient. 

Although most studies of velocity gradients have fo- 
cused on Si II A6355, other features have been examined 
to some degree. While certain features a ppear to have 
little variation from SN to SN (e.g., S II; iPignata et alJ 
|2008[) . Ca II H&K velocity near maximum brightness 
varies significantly amongst SNe la. Additionally, Ca II 
H&K is a strong, blue feature with a large velocity 
gradien t. Using the m odel sp ectra and light curves of 
IKasen fc Plewa (20071) . iFKUl found that the maximum- 
brightness velocity of Ca II H&K strongly correlated with 
maximum-light color. Ca II H&K is an excellent candi- 
date for an alternative way to determine intrinsic color. 



3.1. Sill X6355 



SN la Velocity and Intrinsic Color 



TABLE 2 
Spectral Properties 



SN 



Phase 

(days) 



''Si II 
(10^ kins" 



') 



pEW(Si II) 

(A) 



'JQa II, blue 

(10^ kms-1) 



""Ca II, red 

(10^ kms-1) 



pEW{Ca II) 

(A) 



Ref. 



1981B 
1981B 
1981B 
1981B 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1986G 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 
1989B 



-1.89 

15.00 

17.98 

21.96 

-4.82 

-3.82 

-3.12 

-2.82 

-2.32 

-1.33 

-0.83 

-0.33 

0.17 

0.67 

1.17 

-7.21 

-1.37 

2.60 

4.56 

6.57 

8.57 

10.56 

11.54 

12.52 

13.53 

15.50 

16.47 

17.47 



-12.05 
-10.83 
-10.39 
-10.00 
-11.38 
-11.55 
-11.30 
-11.15 
-11.20 
-10.91 
-10.86 
-10.13 
-10.47 
-10.51 
-10.47 
-11.40 
-10.64 
-10.65 
-10.06 
-10.17 
-10.04 
-10.04 
-9.56 
-9.56 
-9.43 
-9.19 
-9.43 
-9.31 



(0.06) 
(0.06) 
(0.06) 
(0.06) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 
(0.05) 



128.6 (0.1) 
218.2 (0.1) 
169.4 (0.1) 

252.0 (0.1) 

140.1 (0.1) 
131.0 (0.1) 

134.0 (0.1) 

131.2 (0.1) 
126.4 (0.1) 
124.8 (0.1) 
133.2 (0.1) 

122.7 (0.1) 
137.2 (0.1) 
117.7 (0.1) 
131.2 (0.1) 

121.7 (0.1) 
120.6 (0.1) 

119.4 (0.1) 
131.6 (0.1) 

130.8 (0.1) 

136.1 (0.1) 

144.5 (0.1) 
156.1 (0.1) 

130.6 (0.1) 

156.0 (0.1) 

163.1 (0.1) 
163.5 (0.1) 
177.4 (0.1) 



-17.88 (0 

-16.85 (0 

-15.09 (0 

-15.45 (0 

-15.48 (0 

-14.49 (0 

-15.29 (0 

-15.09 (0 

-15.09 (0 

-14.89 (0 

-14.49 (0 

-14.90 (0 



05) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 
05) 

05) 
05) 
05) 



-15.60 
-13.27 
-13.10 

-16.26 
-15.89 

-16.68 
-17.17 
-17.32 

-16.45 
-15.26 
-14.77 
-15.10 
-12.48 
-11.60 
-11.93 
-11.92 
-12.11 
-11.72 
-12.31 
-11.92 
-11.14 

-11.92 
-11.72 
-12.12 



0.06) 
0.06) 
0.06) 

0.05) 
0.05) 

0.05) 
0.05) 
0.05) 

0.05) 
0.05) 
0.05) 
0.05) 
0.05) 
0.05) 
0.50) 
0.50) 
0.50) 
0.50) 
0.05) 
0.05) 
0.05) 

0.05) 
0.05) 
0.05) 



115.8 

79.6 

74.4 

78.7 

100.1 

110.1 

109.2 
92.3 
90.6 

97.1 

115.0 

103.7 

111.7 

126.1 

108.5 

129.1 

99.7 

99.1 

95.4 

84.7 

79.4 

68.8 

73.1 
69.7 
66.9 



0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 

0.1) 
0.1) 
0.1) 

0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 
0.1) 

0.1) 
0.1) 
0.1) 



Ref eren ces. — (1) IWells et al.l jTW^: (2) ICristiani et al.l 1199^'): (3) ILeibundgut et al.l 
figgT); (4) 'Mazzali ct al. (1993); (5) Blondin & Tonrv' (2007); (6) Gomez & Lopez' (1998"); (7) 

[Jefferv et al. (1992); (8) Mazzali et al. (1995); (9) Turatto et al. (1996); (10) Lcibundgut et al. (1993); 
(11) iKirshner et al.. (1993); (12) Patat et aL (1996); (13) SUSPECT, no r eference; (14) Blondin 
et aL, in prep ; (15) Salvo et al. (2001); (16) Li et al. (1999); (17) iTiiFa tto et al. (1998); (18) 

LMatheson et al.i (12008); (19) Branch et al. (2003); (20) Jha et al. (1999F); (21) Blondin et al. (201lJ); 
(22) Garavini et al. (2004); (23) Garavini et al. (2005); (24) Unknown; (25) Garnavich et al. (2001); 
(26) Harnuv et al. (2002); (27) Valentini et al. (2003); (28) Krisciunas et al. (2011); (29) Sauer e t ahl 
(2008); (30) Wang et al. (2003); (31) Benetti et al. (2004); (32) Pignata et al. (2008); (33) Kotak e t^Tl 
(2005); (34) Elias-Rosa et al. (2006b); (35) Stanishev et al. (2007); (36) Anupama et al. (2005); (37) 

[Leloudas et al. (2009); (38) Krisciunas et al. (2007); (39) Altavilla et al. (2007); (40) Pastorello e t~al] 
(2007); (41) Taubenb erger et al. (2008); (42) Wan g et al. (2009b); (43) Garavini et al. (2007b); (44) 

iQuimbv et all (12006^ 1; (451 IKasliwal et al.l (12003) 



In Figure [1] the Si II velocity is shown as a func- 
tion of phase relative to B maximum for the 141 SNe la 
with a velocity measurement, a time of maximum, and 
1 < ATOi5(i?) < 1.5 mag. The SNe with higher velocity 
near maximum brightness also have larger velocity gradi- 
ents. However, at very early times (t < — 10 d), some low 
velocity gradient SNe have similar velocities to high ve- 
locity gradient SNe. Similarly, at later times {t > 15 d), 
it is difficult to distinguish SNe with different velocity 
gradients with a single velocity measurement. Near max- 
imum brightness, the velocity evolution is close to linear 
in time. 

There are two SNe, SNe 2003W and SN 2004dt, that 
are outliers in both early-time velocity and velocity gra- 
dient. SN 2004dt has very high early-time ?;si ij. The 
velocity gradient at early times is consistent with other 
SNe, but also changes dramatically. It declines by 
2750kms~^ in only 1 .4 d, which was previously noted 
(| Altavilla et al.l 120071) . Additionally, the Si II A6355 
feature is wide, flat- bottomed, and potentially con- 
tains two components (jWang et al.ll2006l : lAltavilla et alJ 
I2007f ). SN 2004dt is also a significant outlier in the rela- 
tion between velocity gradient and nebular velocity shifts 
(fMacda ct al...2GlG. ). Finally, SN 2004dt has extremely 



high line polarization ([Leonard et al.l [20051 : iWang et al.l 
l2006f ). In summary, it appears that SN 2004dt is spec- 
troscopically distinct from the majority of SNe la and 
potentially different from even very high velocity SNe. 

SN 2003W has both the highest f si ii measurements of 
any SN in the Fll sample, exceeding —20,000 kms~^, 
and a dramatic velocity decline near maximum light, de- 
clining by over 4000 km s~^ over 5 d (where the typical 
value is 100-200 kms~^ d~^). Its Si II profile appears 
to be flat-bottomed and probably is the result of two 
strong components. Unfortunately, there never appears 
to be two distinct components in any of our spectra. 
Perhaps additiona l spectroscopy at slightl y later times 
would reveal this. iWang fc Wheeler! (|2008[ ) reports that 
there is an unpublished epoch of spectropolarimetry for 
SN 2003W. These data may be extremely informative to 
understanding the nature of this SN, particularly given 
the extraordinary polarization of SN 2004dt. Although 
SN 2003W is clearly a high-velocity and high velocity 
gradient SN, we have excluded it from further analysis 
because of its peculiar velocity evolution. It is unclear if 
SN 2003W is a truly peculiar SN or if the velocity evo- 
lution of other SNe la with extremely high wgi u at early 
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TABLE 3 
Data for Telegraphic Supernovae 



SN 



vgi II 
(kms-i) 



UT Date 



Spec. 
Ref.^' 



MJD 
of tB„.., 



Ref.b 



Phase 



,,0 c 
''Si 11 
(days) (kms~-^) 



1992ag 


-11.8 


19920703.042 


1 


48807.15 




-0.6 


-11.7 


1992al 


-11.7 


19920729.3 


2 


48838.24 




-5.7 


-11.3 


1999cp 


-14.1 


19990619.2 


3 


51363.27 




-15.2 




1999dk 


-16.2 


19990815.06 


4 


51413.33 




-9.4 




2000ca 


-11.0 


20000429.3 


5 


51666.63 




-3.1 


-10.8 


2000dr 


-12.5 


20001008.1 


6 


51833.98 


2 


-8.3 




2001ba 


-10.9 


20010430.17 


7 


52034.36 


1 


-4.7 


-10.6 


2001cj 


-10.5 


20010601 


8 


52066.04 


2 


-4.4 


-10.4 


2001dl 


-10.0 


20010808 


9 


52131.47 


2 


-1.9 


-9.9 


2003gt 


-11.0 


20030805 


10 


52862.17 


2 


-5.6 


-10.7 


2005dc 


-13.3 


20050804.34 


11 


53598.85 


2 


-12.2 




2005ms 


-13.3 


20051231 


12 


53743.45 


1 


—7.7 




2006dm 


-12.0 


20060705.1 


13 


53928.96 


2 


-7.3 




2006ef 


-12.0 


20060824 


14 


53967.95 




1.5 


-12.2 


2006ej 


-13.0 


20060824 


14 


53973.09 




-4.7 


-12.5 


2006en 


-10.2 


20060828.9 


15 


53971.66 




3.7 


-10.4 


2006OS 


-12.7 


20061122.15 


16 


54061.68 




-0.2 


-12.7 


2006qo 


-11.1 


20061201 


17 


54081.96 




-11.1 




2006td 


-10.9 


20061228.77 


18 


54098.92 




-1.4 


-10.8 


2007O 


-10.0 


20070122 


19 


54123.80 




-1.3 


-10.0 


2007sr 


-12.1 


20071220.22 


20 


54448.34 




6.1 


-12.9 


2008dt 


-14.0 


20080701.21 


21 


54646.73 


2 


1.7 


-14.3 


2008ec 


-12.6 


20080716.08 


22 


54674.28 


2 


-10.6 





Spectroscopy references: (1) 'Maza et al.I (11992); (2) 'McNaught et al.' (1993); 
(3}_Jhaetal. (1999a); (4) Salvo et al. (1999); (5) AldcriiiK & Coiilev (2000); (6i 
ISuntzeff fc S mith (2000); (7) Nugent & Wang (2001); (8) Wang (2001); (9) Patat et aQ 
(120011); floT Filippenko et al. (2003); (11) Foley ct al. (2005); (12) Leonard (2005); 
(13 J5eli et al. (2006); (14) Foley et al. (2006); (15) Elias-Rosa et al. (2006a); (16) 
Euhnby et al. (2006a); (17) Silverman et al. (2006); (18) GurugubeUi et al. (2006); 
(191 Silverm an et al. (2007); (20) U mbriaco et al.. (.2007.') ; (211 .Blondin fc Bcrlina 
(|2008l ): (22) Harutvunvan et al. (2008). 

^ Photometry references: (1) Hicken et al.l l l2009al ); (2) IGaneshalineam et aLl (120101 ). 
'^ Uncertainty set to 300 kms~^, which is slightly larger than our expected uncertainty 
for Literature SNe, and larger than the median uncertainty for the Fll sample. 



times is similar to that of SN 2003W. 

A line was fit to the velocity data for each SN (except 
for SNe 2003W and 2004dt) with >3 measurements in 
the window —6 < t < 10 d regardless of A77ii5(i3). We 
fit for 

v{t)^vt + v°, (1) 

where v{t) is the velocity, t is the time relative to B 
maximum, and v and v^ are a slope and offset for the line, 
corresponding to the velocity gradient and velocity at 
time of B maximum (t = 0), respectively. The resulting 
fit values for v and v'^ are shown in Figure [H 

SNe la from the different light-curve shape subsam- 
ples have different relations between Si II A6355 velocity 
and velocity gradient. For SNe la with 1 < Ami5{B) < 
1.5 mag, there is a clear relation between the two val- 
ues such that SNe with a higher velocity at maximum 
also have a higher velocity gradient. This is quantitative 
confirmation of the previously noted correlation. The 
fast decliners tend to have low velocities, but high veloc- 
ity gradients rel ative to the other SN e, confirming previ- 
ous results fe.g-. IBenetti et aLll2005[) . The slow decliners 
have a large scatter, but also have a similar trend to the 
moderate decliners (although they may have a slightly 
lower velocity gradient for a given velocity). 

For the moderate decliners, a line was then fit to the 
values of v and w° , 

v = av" + p. (2) 

Substituting for u, we can construct a family of func- 



tions that describe the velocity evolution near maximum 
brightness and only depend on u*^. 



v{t) 



vt + 


«" 


1 


+ /3)t-f«" 


(1 + 


o£)v^ + ^t 



(3) 



Using the fit values, the family of functions correspond 
to 

i;si II («§; II, i) = 4 ii(l - 0.0322i) - 0.285i, (4) 

where Ug; n is the velocity at maximum, t is measured in 
days relative to maximum brightness, and usi n and Wg; n 
are measured in 10^ km s~^. This family of functions 
describes the velocity evolution of SNe la near maximum 
brightness, and a number of the functions are plotted in 
Figure [Ij 

A test of how well this family of functions performs 
over our chosen time interval (—6 < t < 10 d) is to ex- 
amine the residuals of the functions relative to the data. 
Every velocity measurement within our time interval pro- 
vides an estimate of Wg; n for that SN. For all SNe with a 
velocity measurement within a day of maximum bright- 
ness, we measure the "residual" : the difference between 
the predicted Ugj jj from a single measurement and the 
measured velocity within 1 day of maximum brightness. 
We call this the residual, but it is not exactly the resid- 
ual since the measurements were not obtained exactly at 
i = d. The residuals are shown in the bottom-left panel 
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Fig. 1. — Top left: Temporal evolution of Si II A6355 velocity, lisi n, for the Fll sample of SN la with 1 < Ami5(B) < 1.5 mag. 
Measurements for the same SN are connected by solid lines. The vertical dashed lines mark the time interval over which we fit the velocity 
evolution. The dotted lines represent a subset of the family of functions that describe the velocity evolution (Equation (JJ . Each SN is 
color-coded based on the estimate of iigj jj (using Equation [5j from the measurement closest to t = d. The color bar at the top shows 
how the colors correspond to fSj jj. SNe with no measurement in —6 < t < 10 d are plotted in black. SNe 2003W and 2004dt are plotted 
in black and are excluded from further analysis. Individual measurements from the Telegraphic sample are presented as black points (see 
Section l2.2.3l l. Top right: Histograms of dS- jj. The grey histogram is for SNe where we have measured velocities from spectra. The empty 
histogram represents SNe 2003W and 2004dt. The black histogram represents the SNe from the Telegraphic sample. The horizontal dashed 
line represents the velocity which roughly separated the .W09. sample into "Normal" and "High- Velocity" SNe. Bottom-left: Residuals of 
"si 11 Estimates relative to measured velocity within one day of t = d. The dotted horizontal lines represent the standard deviation of 
220 km s~^ . Bottom-right: Histogram of the residuals. The dotted horizontal line is at 0. 



of Figure [T] 

The residuals are small and centered on zero. The 
average for all residuals is —10 km s^^ and the standard 



deviation is 220 km s~^. There are few residuals that 
are more than 1000 km s~^ from zero. The residuals 
are slightly larger at the extremes of our time interval 
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Fig. 2. — Best-fit values to the slope and offset for a line (Equa- 
tion[TJ describing the vgi n evolution of the Fll sample. SNe with 
Amis < Ij 1 < Amis < 1-5, and Amis > 1-5 mag, are represented 
by blue diamonds, black circles, and red squares, respectively. The 
size of the points is inversely proportional to the size of the un- 
certainty. Each point with 1 < Am,is < 1.5 mag has a circle sur- 
rounding it, with a minimum size to help the reader see the smallest 
points (corresponding to SNe with the largest uncertainty). The 
solid line is the best-fit line for the 1 < Amis < 1-5 mag subsam- 
ple. 



than near maximum brightness, but the difference is still 
relatively small. The typical measurement uncertainty 
is <100 kms^^. Velocity uncertainty from host-galaxy 
rotation is ^300 kms~^. Therefore, the scatter in the 
derived relation is not the dominant source of uncertainty 
for Ugj jj. Combining these sources of uncertainty, the 
total uncertainty is ^400 kms~^. 

Using Equation|4]and a single velocity measurement in 
the time interval —6 < i < 10 d, an estimate of Wg; jj can 
be made for a SN. This allows for a direct comparison of 
all SNe with such measurements. Solving Equation H] for 



V, 





Si II' 



one finds 



^Si II 



= (wsi II + 0.285i)/(l - 0.0322t). 



(5) 



Table [T] lists Ug; jj for each SN in the Fll sample with 
a measurement in —6 < i < 10 d. The value given 
is for the measurement closest to t = d. In Fig- 
ure [31 we plot Wgj JJ as a function of Ami5{B) for the 
full Fll sample. We note that these values are derived 
from Equation [SJ and the velocity estimated for SNe 
with Atoi5(_B) < 1 mag or A?7ii5(i3) > 1.5 mag may 
have systematic offsets (but typically < 1000 kms~^) 
from their true value. Regardless, there appears to be a 
dearth of high- velocity/fast-declining SNe la, but there 
are not many SNe with large Ami^{B) in the Fll sam- 
ple. Slow-declining SNe la can have high-velocity ejecta, 
but there also appear to be a significant population of 
low- velocity/slow-declining SNe la. 

The lack of high- velocity/fast-declining SNe la is not 
unexpected. SNe la with Ami5{B) > 1.5 mag are spec- 
troscopically distinct from slower-declin ing SNe la, being 
spectr oscopically similar to SN 1991bg (|Filippenko et al.l 
Il992at iLeibundgut et al.l \T993.) . No high- velocity 







1.0 1.5 2.0 

Am|5(B) (mag) 

Fig. 3. — Si II velocity at maximum brightness {v^^ jj) as a func- 
tion of Am,is(-B) for the Fll sample. The vertical dotted lines 
mark the region of our full analysis (1 < Ami5(S) < 1.5 mag). 



SN 1991bg-like SN la has yet been discovered. It is worth 
noting that despite having a low velocity near maximum 
brightness, SNe la with Atoir(.B) > 1.5 mag tend to have 
high- velocity gradients (jBenetti et al.l 120051 : Figure [2]) . 
Therefore, Equation |4] may not properly describe the ve- 
locity evolution of these SNe. 

Since fast-declining SNe la have ejecta velocities 
similar to or lower than slower-declining SNe la at 
maximum brightness, one cannot attribute their fast 
post-maxim um decline to a rapidly expanding photo - 
sphere fe.g.. lH6flich et al.lll996HPinto fc Eastmanll2001| ). 
Rather, the post-maximum decline is likely set by the 
amount of Fe-group elements in the ejecta (and thus 
linked to ^^Ni pr oduction in the explosion and the SN 
p eak luminosity ' ) (iKasen fc Woos lcv 2007). 

Be netti et al.l ()2005D also showed a lack of high- velocity 
gradient SNe la with Amic,{B) < 1 mag. Many of 
these SNe are spectros copically similar to SN 1991T 
(jFilippenko et al.lll992 b D. w hich have a low velocity near 
maximum brightness ( W09( ). Although the relation be- 
tween velocity at maximum brightness and velocity gra- 
dient for SN 1991T-like SNe appears to be consistent 
with that found in Equation |4l there are larger residuals 
for these SNe. 

Within the limited range of 1 < Ami5{B) < 1.5 mag, 
there is no clear trend between Wgj jj and ATOi5(i?); how- 
ever, the SNe with the slowest ejecta are found among 
the slowest decliners. 

In Figure JH we show the Wg; jj cumulative distribu- 
tion functions (CDFs) for SNe la in the Fll sample with 
Ami5(B) < 1, 1 < Ami5{B) < 1.5, and Ami5(B) > 
1.5 mag. Using these broad groups, it is easy to see 
the trends described above: there is a lack of high- 
velocity /fast-declining SNe la, there are a significant 
number of high- velocity/slow-declining SNe la, and there 
are a significant number of low-velocity/slow-declining 
SNe la. Additionally, the figure shows that SNe la with 
1 < ATOi5(i?) < 1.5 mag tend to have a higher Ugj jj than 
SNe la in the other groups. Performing Kolmogorov- 
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Fig. 4. — Si II velocity at maximum brightness (tigj jj) CDFs for 
the Fll sample with Ami5(B) < 1, 1 < Ami5(_B) < 1.5, and 
Ami5(_B) > 1.5 mag (dashed blue, solid black, and dotted red 
lines, respectively). 



Smirnov (K-S) tests, we find that it is unlikely that 
SNe la with 1 < /S.mi^{B) < 1.5 mag and those with 
ATOi5(i?) < 1 mag and those with ATOi5(i?) > 1.5 mag 
come from the same parent population {p = 0.0065 and 
0.04, respectively). 

Because of the spectroscopic difference between SNe la 
with 1 < Ar7ii5(i?) < 1.5 mag and many of the SNe la 
outside that range, we find it prudent to restrict our 
analysis to this range. Additionally, the differences in 
velocity populations for the three groups above and the 
different relation between velocity and velocity gradient 
for SN 1991bg-like SNe la (see Figure ^ give further 
reason t o focu s on the limited range in light-curve shape. 
Finally, IW09I found that "Normal" and "High- Velocity" 
SNe la have similar light-curve shape and host-galaxy 
morp hology dis t ributio ns over this range. 

Hic ken et al.l (|2009al ) compiled morphology classifica- 
tions for most of the host galaxies of SNe la in the 
Fll sample. Figure [5] presents v^^ jj as a function of 
host-galaxy morphology for the F ll sample w ith 1 < 
Atoi5(B) < 1.5 mag. As found bv I Wang et "all (|2009al )~ 
but now using full velocity information, we do not see any 
significant differences in the Si 11 velocities observed in 
host galaxies of different morphologies. Although there 
are few very low velocity (wg; u > —10,000 kms~^) 
SNe la in elliptical galaxies, and no high- velocity (wg; jj < 
— 11,800 kms~^) SNe la in galaxies later than Sc, the 
number of SNe in these morphology bins is small. This 
may come from the combination of a slight trend between 
Ug; JJ and Ami5(i?) (Figure |31) as well as faster declining 
SNe la having a h igher frequency in early-type galaxies 
fe.g.. lHowel]|l200l . 

We also examine the Bmax — Knax CDF for sepa- 
rate Wg; JJ bins. Figure [5] shows these CDFs for four 
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Fig. 5. — Si II velocity at maximum brightness (iigj jj) as 
a function of host-galaxy morphology for the Fll sample with 
1 < Am,i5(_B) < 1.5 mag. The black circles represent individ- 
ual SNe, while the red diamonds represent the median value for 
each morphology bin with the error bars representing the median 
absolute deviation for each morphology bin. 

SNe, respectively). There is a general trend for the me- 
dian color to become redder with velocity. While the 
highest and lowest velocity bins have similar shapes and 
appear to simply be offset in color space, the middle two 
velocity bins have different shapes. The second veloc- 
ity bin (-10,800 > ■ug; jj > -11,800 kms'^) has a 
broad distribution of colors, containing the bluest SNe 
and a significant population of highly reddened SNe. 
Meanwhile, the third velocity bin (—11,800 > Wg; jj > 
— 12, 800 kms~^) has a relatively narrow range of colors 
and not as many very red SNe as the second velocity 
bin. Performing K-S tests, we find that the parent pop- 
ulation of the lowest-velocity bin is significantly different 
from the other velocity bins (all with probability <0.04). 
Despite the apparently different shapes, K-S tests do not 
find a statistically significant difference between any of 
the other velocity bins. 

3.2. CallH&K 

In this section, we repeat much of the same analysis 
presented in Section [3TT1 but for the Ca II H&K feature. 
As noted in Section[Tl Ca II H&K is a very strong feature 
that is relatively blue and thus visible in optical spectra 
at z > 1 and shows significant velocity evolution and 
diversity. These qualities make Ca II H&K appealing for 
use as an intrinsic color indicator at high z. 

However, unlike Si II A6355, Ca II H&K is in a com- 
phcated spectral region. Whereas Si II A6355 is a sin- 
gle, isolated feature, other lines are blended with Ca II 
H&K, and the pseudo-continuum near the feature is set 
by a combination of line blanketing and nearby strong 
features. Furthermore, Ca II H&K often shows a compli- 
cated absorption profile, with a flat bottom, shoulders, 
and/or multiple distinct absorption minima, complicat- 
ing any analysis (Figure [7]). 

The velocity of each absorption minimum in the Ca II 
H&K feature of the smoothed spectra is automatically 
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-10.8, -10.8 > vO. jj 



Bmax — Vinax CDFs for the Fll sample with Dgj jj > 
> -11.8, -11.8 > -ug; „ > -12.8, and 
Dg; JJ < -12.8x10'^ kms~^ (dotted black, dot-dashed blue, dashed 
green, and solid red lines, respectively). IW09I set —11,800 kms"^ 
as the dividing line between the two velocity groups outlined in 
that work. 



recorded (Blondin et al., in prep.). We only exam- 
ine spectra with one or two minima. If two minima 
are found, the higher/lower velocities are classified as 
"blue" /"red." If only one minimum is found, it is cate- 
gorized as the red or lower- velocity component. In Fig- 
ure [SI uca H&K found from the blue and red absorption 
minima are shown as a function of vsi ii. Examples can 
be seen in Figure [T) The two absorptions show distinct 
"clouds" representing a clear correlation between wgi ii 
and iJca H&K for each component independently. How- 
ever, there are several red velocity measurements which 
overlap with the blue cloud of points. The most likely 
explanation is that these measurement correspond to a 
higher velocity component that was incorrectly identi- 
fied as a red component. The obvious reason for this is if 
the blue component dominates over the red component 
to the point where there is no distinct minimum for the 
redder component. 

Our goal of connecting velocities and colors will be un- 
dermined if we misidentify the absorption components. 
When usi ii has been measured, there is no ambiguity, 
but for z > 0.4, usi ii is rarely measurable. To remove 
interlopers from our Ca II H&K velocities, we exam- 
ine 4 categories of objects with the properties: /ca = 
fxivca H&K + 9000 km s"^)//A(wca h&k) < 1-5, Ca II 
H&K feature has a visible red shoulder, pEW(Ca II) > 
150 A, and iiCa h&k > -14,000 kms"^ 

A ratio of fluxes in the spectra at vca h&k and offset 
by +9000 kms^^, which we call /ca, is a strong discrim- 
inant. If vca H&K is matched to a blue feature, the flux 
ratio may be small if there is an additional absorption 
component to the red. Conversely, a red feature will gen- 
erally have much higher flux 9000 kms~^ redward of the 
absorption. One of us (RJF) visually inspected (multiple 
inspections were performed for all spectra, with at least 3 
classifications per spectrum, and requiring consistent re- 
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Fig. 7.— Optical spectra of SN 20061e showing the Ca II H&K 
feature in velocity space relative to 3945 A. Their phase is marked 
to the right. The blue and red absorption velocities are shown as 
dashed blue and red vertical lines, respectively. The spectra show 
the range of possible spectral shapes (from top to bottom); a red 
shoulder that the automatic procedure does not detect, causing 
the blue feature to be classified as red, a double-absorption profile 
where the blue component is stronger, a double-absorption profile 
where the two components are of comparable strength, a double- 
absorption profile where the red component is stronger, and a single 
component where no blue absorption is found. 



suits) all spectra to determine qualitative characteristics 
of the spectra. After making all qualitative assessments, 
it was found that those with a noticeable red shoulder 
on the feature were exclusively in the blue cloud (Fig- 
ure [5|) . These spectra correspond to the case where the 
blue component is dominant and the automatic proce- 
dure could not identify the other component. Addition- 
ally, we found that spectra with a large pEW(Ca II) gen- 
erally overlapped with the blue cloud; however, many 
of these spectra were at the highest velocities (for both 
Ca II and Si II) , and their true association was ambigu- 
ous. Since these spectra were also typically at the earliest 
phases (i ^ — 7 d), excluding them will make little differ- 
ence to our final analysis. We therefore chose to include 
these spectra. Finally, very few blue absorption features 
are at wca h&k > 14,000 kms~^. We therefore declare 
any red t^ca h&k with a velocity below this threshold 
to be a true measurement of the red ejecta velocity re- 
gardless of other measurements. Ultimately, spectra with 
vca H&K < —14,000 kms~^ and either if /ca < 1-5 or 
if there was a visible red shoulder, are assumed to be a 
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analysis and are marked with a red 'X.' 



misidentification and were rejected. Graphically, these 
samples can be seen in Figure [H 

In Figure [9l the red Ca II velocity is shown as a func- 
tion of phase relative to B maximum for the 685 spectra 
of 109 SNe la with a Ca II H&K measurement, a time 
of maximum, and 1 < Ami^{B) < 1.5 mag. Only spec- 
tra which have passed our previous cuts are shown. The 
velocity evolution of Ca II is similar to that of Si II, but 
Ca II is typically at a higher velocity at a given epoch 
(see Figure [SJ. 

There are fewer spectra and fewer SNe with a wca h&k 
measurement than for wsi ii- This is because (1) the CfA 
FAST spectra extend only to ^3500 A in the blue, and 
the Ca II H&K absorption minimum is often blueward 
of these spectra, and (2) the above procedure rejected 
many spectra. Regardless, the same trend that higher 
velocity SNe have higher velocity gradients is clear. This 
is also shown in Figure I10| which is similar to Figure [2l 

It appears that early-time {t < —6 d) Ca II H&K data 
might be a cleaner diagnostic of maximum-light velocity 
and velocity gradient than for Si II A6355, but the small 
number of spectra at these phases prevents a definitive 
statement. This might be the result of inspecting and 
rejecting individual spectra for Ca II H&K and not for 
Si II A6355, which might also have multiple components 
at early phases. SNe with high-velocity Ca II H&K fea- 
tures continue to have higher velocity than their lower 
velocity counterparts to later phases than is seen with 
Si II A6355. Despite this, the linear velocity evolution 
does not maintain at these later phases (there appears 



to be a plateau for most SNe). As a result of a paucity 
of early-time data and the changing velocity gradient at 
later times, we define the phase range of— 4<t<9d 
to be the range over which Vq^ h&k '^^^ ^^ cleanly mea- 
sured. 

Performing the same procedure as above (again, ig- 
noring SN 2004dt, which also has peculiar Ca II H&K 
evolution — SN 2003W only has a Ca II H&K velocity 
for our i = — 11 d spectrum, so it is also not included) 
using uca H&K for SNe la with 1 < A?7ii5(_B) < 1.5 mag, 
we determine the relation 

WCa H&K (Wca H&K>0 = ^Ca H&k(1 - 0.0515^) - 0.633i. 

(6) 

Examples of this family of functions are shown in Fig- 
ure [HI 

The residuals (as defined in Section lOT) . have a larger 
scatter {a = 570 kms^^) than that of Si II A6355, but 
the residuals are not that large considering other fac- 
tors contributing to the total uncertainty in this mea- 
surement. It is unclear if the larger scatter is the result 
of a poorer representation of the velocity gradient be- 
cause of less data, less linear velocity evolution of Ca II 
H&K, a less linear relation between velocity at maxi- 
mum brightness and velocity gradient, interloper mea- 
surements, and/or the relation has larger intrinsic scat- 
ter. There are also more "catastrophic" outliers (resid- 
uals more than 1000 kms~^), which are likely the re- 
maining absorption velocities that overlap with the blue 
cloud. 

Similar to what was done with Equation |4] and Si II 
Ca H&K '^^^ be found using Equation [6l 



A6355, V? 



"Ca H&K 



(«Ca 



H&K + 0.631t)/(l-0.0513t). (7) 



Table [T] lists v^^ jj^j^ for each SN in the Fll sample 
with a measurement in —4 < t < 12 d. The value given 
is for the measurement closest to i = d. 

4. PSEUDO-EQUIVALENT WIDTH 

Although IFKIII showed that ejecta velocity and in- 
trinsic color were connected and we show in Section [5] 
that the relation is robust, ejecta vel ocity m ay not be 
the best indicator of intrinsic color. iFKllI postulated 
that the redder colors in high-velocity SNe may be the 
result of broader lines causing increased line blanketing 
in the near-UV. If this explanation is correct, one might 
expect line widths or line strength to be more strongly 
correlated with intrinsic color. 

FWHM and pEW measurements have observational 
advantages over velocity measurements. Specifically, 
even if the redshift of a SN la is poorly known (e.g., 
no host galaxy redshift for high-z SNe), the FWHM and 
pEW is still robust. Unfortunate ly, pEW is sensitiv e to 
galaxy light contamination (e.g.. iFolev et al.ll2008al) . If 
FWHM is measured as part of fitting a Gaussian profile 
to a spectral feature, then deviations from that profile 
will cause the measured FWHM to lack a physical con- 
nection to the data. 

For the Fll sample, the measured Si II A6355 FWHMs 
assume a Gaussian profile. This has not been attempted 
for the Ca II H&K feature since it rarely has a Gaus- 
sian profile. We performed the below analysis for the 
FWHM(Si II) measurements, but found pEW(Si II) mea- 
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Fig. 9. — Same as Figure [T] except for iiCa 



H&K- 



surements to yield a more robust result. We therefore (|2005f ) classifications based on velocity gradient 

only report the pEW results below. 

Several groups have examined pEWs of Ca II HfcK 
and Si II A6355 at bot h low a n d high z (e.g., iFolatellil 
200l iHachinger et all [200l iGaravini et all 12007a : 



4.1. 



Broiider et al.ri2008t iFolev et all l2008at iBalland et all 



200^ 'Branch et al. ■2009; 'Blondin et al .l 1201 Ibl: 



Konis hi et al. 2011; Nordin et al. 2011a; W alker et aP 
20Tl . IBranch et all (|2006h used the pEW of Si II A6355 



along with the pEW of Si II A5972 to subcl assify SNe la 
His classifications match up well with the iBenetti et al.l 



Si II A6355 Pseudo-Equivalent Width 

As noted above, several studies have examined pEWs 
for several features, including Si II. In Figure 111) we 
show how pEW(Si II) is related to the spectral phase, 
Amic,(B), and vsi ii. There is a wealth of information in 
these data, but we will only briefly discuss them, leaving 
further analysis to future work. 

Examining Figure lll[ we note several observational 
facts. First, as expected, higher-velocity SNe tend to 
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Fig. II. — pEW(Si II) as a function of vgi n for different pliase 
and Am,i5{B) ranges. Tlie panels (clockwise from lower-left) show 
a phase range of 3 < i < 10, —10 < t < —3, —3 < t < 3, and 
— 10 < t < 10 d, respectively. The blue diamonds and red squares 
represent SNe la with Amir^{B) < 1 and Ami5{B) > 1.5 mag, 
respectively. The dark grey, black, and light grey circles repre- 
sent SNe la with 1 < Ami5(S) < 1.5 mag and -10 < t < -3, 
—3 < t < 3, and 3 < t < 10 d, respectively. The size of the sym- 
bol corresponds to the inverse of the uncertainty. A solid line is 
plotted to correspond roughly to the relation between vgi n and 
pEW(Si II) at early times and high velocity. The same line is 
plotted in the upper-right and lower-left panels to aid the eye in 
comparing the different phases. The lower-right panel shows the 
moderate decliners (1 < Ami5(B) < 1.5 mag) for the three phase 
ranges. 



have stronger lines. Second, slower-declining SNe tend 
to have weaker lines. Third, lines tend to get stronger 
over the phase range shown here (—10 < t < 10 d). 
We also note a generic trend where SNe la with t^si n > 
— 12,000 kms~^ can have a large range of pEW that 
appears to be completely unrelated to the measured ve- 



locity. However, above this velocity, there appears to be 
a strong relation between pEW(Si II) and wsi ii at all 
phases shown here. This is likely the result of the line 
saturating. This effect may also be connected to the fact 
that SNe la with wg; jj > -12,000 kms"^ have a small 
scatter in intrinsic color, while high er velo city SNe la 
have a larger range in intrinsic color (jFKllt Section [5]) . 
Figure [12] shows the pEW evolution as a function of 
phase for individual SNe. Similar to the method outlined 
in Section [31 we find a family of functions which describe 
the pEW evolution as a function of time. Using the same 
process as before, we can determine pEW(i = 0) = pEWo 
from a single pEW measurement in the appropriate phase 
range (-7 < t < 7 d for pEW(Si II)). We find 

pEWo(Si II) = (pEW(Si II) + 4.48t)/(l - 0.0339i). (8) 

In the phase range — 7 < i < 7 d, the evolution is rel- 
atively linear. For most SNe, the pEW changes <20% 
over this phase range; therefore, knowing the exact evo- 
lution of this feature with time is not critical. Unlike 
T;g; JJ, pEWo(Si II) has a relatively Gaussian distribution 
for our sample. 

Table[I]lists pEWo(Si II) for each SN in the Fll sample 
with a measurement in— 7<i<7d. The value given is 
for the measurement closest to i = d. 

4.2. Ca II H&K Pseudo-Equivalent Width 

Unlike wca h&k, pEW(Ca II) is not systematically af- 
fected by multiple absorption components. Although in- 
correct definition of the borders of the feature can pro- 
duce incorrect measurements, misidentifying a blue ab- 
sorption component as a red absorption component does 
not change any measurements. 

Figure [T3| shows the Ca II H&K pEW evolution as a 
function of phase for individual SNe. Using the same 
method outlined in Section 14. 1[ we fit a family of func- 
tions to the near-maximum (— 7 < i < 9 d) temporal evo- 
lution of pEW(Ca II), which provides a way to convert 
pEW measurements in that phase range to a maximum- 
light measurement, pEWo- We find 

pEWo(Ca H) = (pEW(Ca H) + 0.708t)/(l - 0.0210t). 

(9) 
Similar to Si II A6355, Ca II H&K has a relatively 
Gaussian pEWo distribution for our sample. 

Table [U lists pEWo(Ca II) for each SN in the Fll sam- 
ple with a measurement in— 7<i<9d. The value 
given is for the measurement closest to t — Q d. 

5. INTRINSIC COLOR 
5.1. Matching the Fll Sample to the W09 Sample 

The 'WOg' SN la sample is large. It is further useful 
for various studies since it is separated into two groups 
by ejecta velocity. However, it lacks published velocity 
measurements. We are able to provid e velo city measure- 
ments for 59 of the 121 S Ne la in the lW09l sample (after 
applying the same cuts as lFKllt ). 

WP9 used the Si II A6355 velocity evolution of 10 well- 
observed low-velocity gradient SNe la to determine an 
average velocity as a function of phase and a disper- 
sion around that average. SNe which had velocity mea- 
surements >3cr above the average in the time interval 
— 7 < i < 7 d were considered "High- Velocity" SNe, while 
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Fig. 12. — Same as Figure[l] except for pEW(Si II). 



all others with measurements within that time interval 
were considered "Normal." Since there are intermediate- 
velocity SNe and since SN velocities at later times are 
similar, this method produces a coarse assessment of the 
SNe with potential ambiguity. 
Figure [T4l shows histograms of Ug; tt for the Fll sam 



pie of SNe la that overlap with the IW09I sample. The 
two histograms represent SNe classified as "Normal" and 
"High- Velocity" by IW09I . The two subsamples are rel- 
atively separate, but there is a significant amount of 
overlap for -11,000 > wg; „ > -12,000 kms^^. The 



standard deviation for our relation that measures Vg; jj 
is 220 km s^"'^, the typical measurement uncertainty is 
<100 kms^^. (Velocity uncertainty from host-galaxy ro- 
tation does not affect this direct comparison.) Given the 
range of overlap of the subsamples, the data sets and/or 
methodology must be responsible for some of the over- 
lap. This gives additional motivation for using individual 
velocity measurements to describe the ejecta velocity of 
a SN la. 

5.2. Intrinsic Color for the WOO Sample 
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In Figure [151 we plot the light-curve (but not host- 
galaxy redde ning) corrected peak absolute magnitude 
of the IW09I matched sample as a function of their 
Bma x — Kn ax pseudo-color. This is similar to Figure 4 
from lFKTll . The symbols representing each SN are color- 
coded by ■yg; jj. As noted bv lFKlll the higher- velocity 
SNe tend to be redder (brighter) for a given peak mag- 
nitude (color). 



Using Equation 1 from lFKlll . we have a relation be- 
tween measured E{B — V) and the peak magnitude. The 
IW09l E{B - V) and the Fll S^ax - Kiax values are 



highly correlated. Fitting a line to the data, we find a 
simple offset of Bmax ~ Knax — E[B ~V) ~ 0.081 mag. 
Th is is consistent with t he offset (-0.070 ±0.012) found 
bv iPhillips et al.l (|1999( ). especially considering the Fll 
-Bmax — Knax valucs iucludc a Milky Way correction. Ap- 
plying this offset to the best-fit values of Equation 1 of 
IFKIII for the "Normal" SNe, we obtain the solid line in 
Figure [TSl This line represents the dust reddening for a 
SN la with an "intrinsic" color of zero. Although this 
may not be the true intrinsic color, deviations from this 
line in the color direction represent deviations from a 
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Fig. 14. — Histograms of ijgj jj measured from the Fll spectra. 
The blue and red histograms represe nt SN e classified as "Normal" 
and "High Velocity," respecti vely, by lWOSI . The dashed line repre- 
sents the value of Vgj jj that IW09I nominally used to separate the 
two subsamples. 
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Fig. 15. — The light-curve shape corrected peak absolute V 
brightness as a function of Bmax — Vmax. The color of the sym- 
bol corresponds to the SN's Vgj jj, with the color bar at the top 
of the figure displaying the correspondence. The solid black line 
represents the host-galaxy dust reddening of a zero-color SN la. 
Deviations from this line in color is the intrinsic color of the SN. 



nominal color that is corrected for host-galaxy redden- 
ing. Specifically, we define 






where M, 
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M,p)/Rv^C, (10) 
a(Ami5(B) - 1.1), Af,p = 



max ^ max 
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Fig. 16. — Distribution of intrinsic Bmax- 
matched sample of SNe la. 



■Vm 



0.3 

for the lWOSl -Fll 



-19.26 mag, a = 0.75, Ry = 2.5, and C = 0.081 mag 
(iFKllD . This measurement may be offset from the true 
intrinsic color, but the same offset should apply to all 
SNc, making the measurement useful for directly com- 
paring SNe. 

We restrict the following analysis to SNe la with 
Smax — V^nax < 0.4 — 0.081 = 0.319 mag, correspond- 
ing to the eIb -V) < 0.4 mag cut used bv lFKlll . This 
is the color range where low and high- velocit y SNe la 
appear to have the same reddening law (|FKlll ). 

Figure [16] presents a histogram of the intrinsic Bmax — 
Knax for the sample. The shape is slightly skewed to red 
colors, with a skewness of 0.48. The shape of the color 
distribution is similar to the color distri bution expected 
from t he asymmetric explosion models of lKasen fc Plewal 
()200 7ll when accounting for the frequency with which a 
particular viewing angle is observed. 

SN la color is correlated w ith light-curve shape (e.g., 
iRiess et al.lll996l: lTripdll998D. with faster dechning SNe 
having redder colors. iMaeda et all ()2011|) corrected the 



observed colors of 
lation between Bmax 
|R)latelh et al. (2010) 



SNe 



m their sample using the re- 
— Vinax and Ami5{B) found by 

^ ^ Over the full Am^^( B) range used 

for this analysis, the iFolatelli et al.l (|2010f ) relation pre- 
dicts a correction of 0.060 ± 0.025 mag. Although there 
is only a slight correlation between light-curve shape and 
ejecta velocity for the range of A77ii5(_B) examined here 
(e.g.. Figure [3]), there may still be a significant correla- 
tion between Atoi5(_B) and intrinsic color for the IW09l - 
Fll matched sample. Figure [TTl sho ws the distribution 
of color and light-curve shape for the lW09l -Fll matched 
sample. The two quantities displayed (Bmax — Knax and 
Ami5{B)) are not correlated for the sample, having a 
correlation coefficient of 0.04. Therefore, the intrinsic 
color is not significantly affected by light- curve shape for 
the sa mple. Given the uncertainty in the IFolatelli et al.l 
()2010D relation, the fact that it was derived using a com- 
bination of low and high-velocity SNe la, the small cor- 
rection expected for SNe la in this analysis, and the lack 
of a correlation between light-curve shape and intrinsic 
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Fig. 17. — The int rinsic Bmax — Vmax pseudo-color as a function 
of Ami$(B) for the lW09l -Fll sample. The correlation coefficient 
is 0.04. 



color for the sample, we do not apply this correction. 

In Figure [IHl we compare the intrinsic color of each SN 
to its Wg; jj. There is a trend where SNe with redder in- 
trinsic colors systematically have higher ejecta velocity. 
The correlation between the two quantities is 0.28, and a 
linear least-squares fit results in a slope that is significant 
at the 3.4-tT level. Performin g a Bavesian Monte-Carlo 
linear regression on the data ()Kellvll2007[ ). we find that 
99.142% of the realizations have a negative slope. The 
trend is more apparent when separating the sample into 
equal-numbered groups based on Wgj jj and examining 
their median values. The tre nd is also qualitatively sim- 
ilar to that found bvlFKllI for the synthetic spectra of 
IKasen fc Plewal (1200711 

From the IW09I sample, it is unclear if the trend is 
linear or if there is a simple color offset. There is a 
larger intrinsic color scatter for the higher velocity SNe 
than for the lower velocity SNe. Splitting the sample by 
Ug; JJ = —11,800 kms~^, the intrinsic color scatter is 
0.070 and 0.047 mag for the high and low-velo city su b- 
samples, respectively. This was postulated by iFKllI as 
the reason that the "Normal" SNe la produced a smaller 
scatter in Hubble residuals than the "High- Velocity" 
SNe la. However, even after removing a linear trend 
between velocity and intrinsic color, the higher-velocity 
SNe la continue to have a higher intrinsic color scatter 
(0.071 vs. 0.046 mag). This suggests that SNe la with 
lower ejecta velocity may be intrinsically better distance 
indicators, and with enough SNe la, it may be prudent 
to reject high- velocity SNe la from cosmological samples. 

5.3. Intrinsic Color for the Fll Sample 

As describ ed in Section 12.11 there are some drawbacks 
to using the IW09I data. Here we repeat the above anal- 
ysis to determine the intrinsic color of SNe la, but use 
the Fll sample. Specifically, deri ved photometric val- 
ues have be en compiled from CfA3 (iHicken et al.ll2009al) 
and LOSS (jGaneshalingam et al.|[2010D . The CfA3 peak 
magnitudes were corrected to match the LOSS sys- 



tem. See Section 12.11 for details. We derive intrinsic 
-Bmax— Knax pscudo-colors for cach SN in the Fll sample. 
Using the measurements of wsi iij ^'Ca h&k, pEW(Si II), 
and pEW(Ca II) from CfA and literature spectra near 
maximum brightness and the time evolution of these 
quantities as defined by Equations El [71 [HI and [9l we 
derive ^ jj, i;0^h&k; pEWo(Si II), and pEWo(Ca II), 
respectively, for individual SNe. Below, we will compare 
these values to their intrinsic colors. The Fll sample 
takes no data directly from W09i . but the samples have 
some SNe in common and may share derived light-curve 
values. Furthermore, the Wg; jj values used in the above 
analysis are identical to those presented here for the SNe 
in both samples. 

Similar to what was done in Section 15.21 we restrict 
the following analysis to SNe la with iJmax — Knax < 
0.319 mag, the color range where low and high- velocit y 
SNe la appear to have the same reddening law (|FK11D . 
We exclude SNe la with z < 0.01, except f or those with 
Cepheid distances from iRiess et all (|2011l ). There are 
65, 59, 42, and 49 SNe la with measured maximum- 
brightness light-curve parameters, 1 < ATOi5(i?) < 
1.5 mag, and a maximum-brightness spectral parame- 
ter for ^ jj, pEWo(Si II), w[^,h&K' and pEWo(Ca II), 
respectively. 

Using the Wg; jj data. Figure [HI shows the light-curve 
(but not host-galaxy reddening) corrected peak abso- 
lute magnitude of the Fll sample as a function of their 
^max — ^nax pscudo-c olor. This is similar to Figure [T5] 

As was seen in the IW09I sample, the higher velocity 
SNe la in the Fll sample tend to be redder than lower 
velocity SNe la. Using the procedure outlined above, we 
again calculate the intrinsic color for the sample. Fig- 
ure [20l compares the intrinsic color of each SN to its 
ejecta velocity (as derived from Si II A6355 and Ca II 
H&K) and pE\ y valu es. 

As with the IW09I sample, the Fll sample shows a 
strong trend between Ugj jj and intrinsic color such that 
intrinsically redder SNe la tend to have higher velocity 
ejecta. Unlike the W09 sample, the Fll sample has an 
obvious linear relation between these quantities — not 
just an offset, corresponding to 

(S„,ax - Knax)o = (-0.39 ± 0.04) - (0.033 ± 0.004) 

X {vl; ii/lOOO km s"^) mag. (11) 

Performing a linear least-squares fit results in a non-zero 
slope that is significant at the 8.8-cr level. The corre- 
lation between the two quantities is —0.39. Performing 
a Bayesian anal ysis o f Monte-Carlo linear regressions on 
the data (,Kellvll2007f) . we find that 99.949% of the real- 
izations have a negative slope and the median correlation 
coefficient for the realizations is —0.47. There is a rea- 
sonable amount of scatter to the relation, but on average, 
the trend is quite clear. Although slightly offset, the lin- 
ear relati on is s trikingly similar to the theoretical trend 
found bv lFKUl . 

Similar to what was found with the IW09I sample, the 
intrinsic color scatter is larger for higher-velocity (as de- 
termined by Wgj jj) SNe la. The intrinsic color scatter 
for SNe la with t^gj jj higher/lower than —11, 800 kms""'^ 
is 0.095 and 0.072 mag, respectively. After correcting 
for the linear trend, the scatter does not change much, 
resulting in a scatter of 0.093 and 0.070 mag for the sub- 
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Fig. 18. — Maximum-light Si II A6355 velocity (v^^ jj) as a function of intrinsic -Bmax — Vmax pseudo-color for the'W03 sample of SNe la. 
The red circles represent the median values of equal-numbere d velocity bins w i th the error bars being the median absolute deviation. The 
dashed line is the relation between velocity and color for the IKasen &: Plewal II2007I) models (Equation 3 of IFKllT) . The solid line is the 
best-fit linear model for the data. 
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Fig. 19. — Same as Figure [151 except for the Fll sample. 



samples, respectively. However, correcting for the lin- 
ear trend reduces the overall scatter from a = 0.087 to 
0.080 mag. 



Comparing pEWo(Si II) to the intrinsic i?max — V^iax 
pseudo-color (Figure [20]). there is a reasonable correla- 
tion {p = 0.28), with a linear relation defined by 

(S,nax - Kiax)0 = ("0.16 ± 0.03) - (0.0016 ± 0.0003) 
X pEWo(Si II) mag. (12) 

The slope of the best-fit line for the two quantities is 
6.3-cr from zero. Again, looking at medians (now for 
bins of equal number in pEW space), the trend is clear. 
Performing a Bayesia n analysis of Monte-Carlo linear re- 
gressions on the data (jKelly"2007), we find that 98.6% of 
the realizations have a positive slope and a median cor- 
relation coefficient of 0.34. The intrinsic color of SNe la 
can be derived from both the velocity and pEW of Si II 
A6355 near maximum brightness. 
The Ca II H&K measurements produce less robust re- 



sults. Intrinsic color is correlated with 





Ca H&K' 



where 



(similar to Si II A6355) SNe la with higher velocity ejecta 
tend to be intrinsically redder. A linear least-squares fit 
to the data results in a slope that is 4.7-cr from zero. 
However, the data have a large scatter and the correla- 
tion is not apparent by eye (p = —0.24). A Bayesian 
Monte-Carlo analysis of the data results in 96.1% of the 
realizations having a negative slope and a median cor- 
relation coefficient of —0.33. Nonetheless, lower-velocity 
SNe la, as determined from Ca II H&K, appear to have 
less intrinsic color scatter than higher-velocity SNe la. 
Splitting the sample at v^^ ^^-^ = —14,000 kms~^, we 
find that the higher/lower velocity SNe la have intrinsic 
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Fig. 20.— Maximum-light Si II A6355 velocity, pEW(Si II), 
Ca II H&K velocity, and pEW{Ca II) as a function of intrinsic 
Smax — Vmax pscudo-color (top to bottom panel, respectively) for 
the Fll sample of SNe la. The top panel is similar to Figure [TSl 
except for the Fll sample. The red circles represent the median 
values of equal-numbered velocity or pEW bins with the error bars 
representing the median absolute deviation. The dashed lines in 
the to p and third panels are the relation between velocity and color 
for the lKasen FP lcwa (2007) models (Equations 3 and 4 of FKil 
respectively). The theoretical relations are strikingly similar and 
particularly poor for the the Si II A6355 and Ca II H&K data, 
respectively. The solid lines are the best-fit linear model for the 
data. 



color scatter of 0.095 and 0.058 mag, respectively. There- 
fore, a velocity cut may improve cosmological results. 
The linear fit to the data i s significantly d i fferen t from 
that found by|FKll|for the lKasen fcPlewal (pOOTi ) mod- 
els, and detrending the data by the best-fit line slightly 
increases ('-^^0.003 mag) the scatter of the velocity sub- 
samples (but decreases the scatter of the full sample by 
~0.002 mag). 

The pEW(Ca II) measurements do not show a corre- 
lation with intrinsic color. The correlation coefficient is 
0.03, and there is no significant linear trend as found 
by a fit to the data (performing the similar linear fits 
provided a non-zero slope at the 1.8-cr level and positive 
slopes for 65.1% of the realizations). There is no clear 
way to make a cut on pEW(Ca II) to reduce the intrinsic 
color scatter. Two competing effects drive the value of 
pEW(Ca II). First, higher velocity ejecta correspond to 
broader lines (as clearly seen for Si II A6355; Figure [TT|) . 
so SNe la with higher ejecta velocity, which have red- 
der intrinsic colors, should have higher pEW. However, 
SNe la with redder intrinsic B ~ V colors, which have 
higher ejecta velocity, will have a depressed UV contin- 
uum, and should have a lower pEW. These competing 
effects may reduce any correlation for pEW(Ca II). 

6. DISCUSSION & CONCLUSIONS 

SN la kinematics (as measured by the minimum of the 
absorption or the width of the absorption of spectral fea- 
tures) are related to the intrinsic color of a SN la. The 
correlation between SN la ejecta velocity and intrinsic 
color found by FKll' has been further investigated here 
by examining additional kinematic probes: velocity and 
pEW measurements for both Si II A6355 and Ca II H&K. 
Using spectral series of many SNe la, families of func- 
tions, which provide a measurement of velocity and pEW 
at maximum brightness given a measurement near max- 
imum brightness, were constructed. Both features show 
a strong correlation between their velocity gradient and 
velocity at maximum brightness such that SNe la with 
larger velocity gradients tend to have higher velocities at 
maximum brig htness. Similar to high velocity-gradient 
SNe (THaching eTet al.l I2OO60 . higher-velocity SNe also 
tend to have broader features (see also IW09D . Restrict- 
ing the sample to 1 < /S.mi<^[B) < 1.5 mag, we find no 
correlation between host- galaxy morphology and Wg; jj. 
Using the families of functions, SNe la with spectra at 
different phases can be directly compared. These func- 
tions provide a way to directly compare SNe la using a 
single near-maximum spectrum, providing an easy way 
to categorize SNe la in a manner similar to how Attiis 
has been used for photometr y. Using the photom etric 
parameters of |W 09 as well a s iHicken et al.l ()2009aD and 
IGaneshalingam et al.l (|2010[ ) along with the reddening 
law derived in lFKlll we a re able to determine the intrin- 
sic color of both the lW09l and Fll samples, respectively. 

Comparing the maximum-light velocity and pEW mea- 
surements with the intrinsic color data for individual 
SNe, we find correlations between intrinsic color and 



,,0 



,,0 



'^Si II' '^Ca H&K' ^^"^ pEWo(Si II), although the corre- 
lation with v^^ n&iK i^ ^^^^ robust. We find no correla- 
tion between intrinsic color and pEWo(Ca II), and spec- 
ulate that the competing effects of higher- velocity SNe la 
having broader lines and depressed UV continua balance 
each other to some extent, negating any correlation. The 
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intrinsic color of a SN la varies linearly with both Vg^ jj 
and pEWo(Si II), although there is significant scatter. 

Thus far, we have only discussed certain aspects 
of how our measured parameters correlate with each 
other. Specifically, we find that vsi ii and wca h&k 
are highly correlated (Figure [8|). We also find that 
usi II and pEW(Si II) are highly correlated for the high- 
est velocities, but uncorrelated at lower velocities (Fig- 
ure HH). In Figure [H v^, „, pEWo(Si II), t;°^ h&k. 
and pEWo(Ca II) are plotted against each other for 
the Fll sample. SNe with intrinsic color information 
are color-coded by their intrinsic color. Most parame- 
ters are highly correlated with each other (the absolute 
Pearson correlation coefficient for each combination is 
presented in Figure [2T|). The exceptions are Wg; jj and 
pEWo(Ca II), and pEWo(Si II) and pEWo(Ca II). This 
is not unexpected since Wgj jj and pEWo(Si II) correlate 
with intrinsic color while pEWo(Ca II) does not. 

There are a handful of outliers for each relation be- 
tween the spectral parameters, but the outliers do not ap- 
pear to have a strong relation with intrinsic color. There- 
fore, measuring multiple parameters may not directly im- 
prove our knowledge of intrinsic color. Moreover, it is not 
clear if rejecting SNe based on their placement in this 
multidimensional parameter space would improve any of 
the previously determined relations. 

Explosion models combined with radiative transfer 
models must account for the relation between ejecta kine- 
matics and intrinsic color. Accurate SN la models should 
reproduce both the general trends seen in our data as well 
as the scatter. We find that high-velocity (as defined 
by either Wg; jj or v^^ h&k) ^^^ ^^ have a larger scat- 
ter in their intrinsic color even after accounting for the 
velocity-color relation. FKll suggested that the larger 
Hubble scatter seen for "High- Velocity" SNe la is caused 
primarily by the larger scatter in intrinsic color, but also 
suggested that this is primarily driven by the larger ve- 
locity range for those SNe. This appears to be correct, 
but the residual scatter after correcting for velocity indi- 
cates lower-velocity SNe la are better standard crayons 
than higher- velocity SNe la. 

To reduce the scatter and potential bias ([FKllI ) of their 
samples, future cosmological analyses with SNe la may 
elect to reject higher velocity SNe la, along with or in- 
stead of fitting a relation between velocity and color. Fur- 
thermore, the most precise a nd accu rate SN la distances 
must account for this effect. iFKllI outlined how ignor- 
ing this effect will increase the scatter associated with 
SN la distances and may bias those distances, resulting 
in biased measurements of cosmological parameters. 

Although we provide a prescription for determining the 
intrinsic color of SNe la with a rest-frame optical spec- 
trum near maximum brightness, this relation is not ro- 
bust for high-velocity SNe la if the spectrum does not 
cover Si II A6355. For SNe la at z > 0.4, this feature is 
rarely observed. Alternatively, a relatively homogeneous 
data set can be obtained by making a cut based on Ca II 
H&K velocity. SNe la with v^^ h&k ^ -14,000 kms"! 
have a relatively small range of intrinsic Smax — ^ax 
pseudo-color {a — 0.06 mag). To avoid biases and to pro- 
duce the most precise and accurate SN la distances, all 
future SN la cosmology surveys should strongly consider 
obtaining spectroscopy near maximum light for their pri- 



mary cosmological sample. Measuring Si II A6355 ap- 
pears to have some advantages over the easier to ob- 
serve Ca II H&K. Therefore, surveys targeting SNe la at 
z > 0.4 may want to explore NIR spectroscopy. Even 
Ca II H&K is redshifted out of the optical window for 
z > 1.2, making NIR spectroscopy necessary for the 
highest-redshift SNe la. We note that if one has a detec- 
tor that is sensitive to 2 fim, one can in principle measure 
Si II A6355 and Ca II H&K to z « 2.2 and z w 4, respec- 
tively. 

We could not have performed the analysis presented 
here without the large CfA data set. Although fur- 
ther analysis with an expanded sample is necessary to 
strengthen the claims presented here and examine the po- 
tential for combining different measurements to improve 
our knowledge of the intrinsic color, the current sample 
has provided insight and direction for these future stud- 
ies. To expand the current sample, our group will soon 
release its CfA4 sample of SN la light curves. On a sim- 
ilar time scale, the Carnegie Supernova Project (CSP) 
will provide its next sample of light curves. For spec- 
troscopy, the Berkeley Supernova la Program (BSNIP) 
will soon release its s ample of SN la spectra, including 
all spectra used in the lW09l analvsis. Similarly, the CSP 
will soon release their spectroscopic data set. While the 
BSNIP sample will not have the well-sampled spectral 
series of the CfA and CSP samples, it does have spectra 
of many SNe la near maximum brightness, and its spec- 
tra typically cover bluer and redder wavelengths than 
the CfA spectra. Therefore, it should provide many ad- 
ditional SNe for this particular type of analysis. On a 
slightly longer ti me scale, the Palomar Tr ansient Factory 
(iRau et al.ll2009[ ) and Supernova Factory (lAldering et al.l 
|2002| ) will provide light curves and spectra of thousands 
of SNe la. With these combined data sets, it should be 
possible to further refine and extend the results presented 
here. 
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Fig. 21.— Maximum-light Si II A6355 velocity, pEW(Si II), Ca II H&K velocity, and pEW(Ca II) for the Fll sample. SNe with no 
intrinsic color information are black, while those with intrinsic color information are color-coded by their intrinsic color with the mapping 
represented by the color bar. The Pearson correlation coefficient for each set of parameters is noted in each panel. 
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